Excitation Spectra of Correlated Lattice Bosons in a Confining Trap 
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We consider ultracold bosonic atoms in an optical lattice with an external trapping potential. 
To study the excitation energies of the resulting Bose-Hubbard model, we develop a method based 
on a time-dependent generalization of the Gutzwiller ansatz. We calculate the excitations of the 
homogeneous system both in insulating and superfluid regime, concentrating particularly on those 
near the superfluid-Mott insulator boundary. Low-lying excitation energies in presence of a static 
harmonic trap are obtained using this method and compared with the homogeneous case. Further, 
we explore the dynamics of the center of mass and the breathing mode in response to time-dependent 
perturbations of the trap. 
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The study of ultracold bosonic atoms_ confined in an 
optical lattice as a correlated system [H has been an 
active area of research in the past few years 0. It 
could provide a better understanding of many phenom- 
ena driven by strong many-body correlations in a con- 
trolled fashion, which is rather difficult to explore in con- 
ventional solid state systems. These include the quantum 
phase transition between the Mott insulating (MI) and 
the superfluid (SF) phases [3|, |j], Luttinger liquid behav- 
ior in quasi-one dimensional condensates in the Tonks- 
Girardeau regime [5|, etc. It has been proposed that 
these systems could be used to simulate various quan- 
tum spin models jf|, especially in reduced dimensional- 
ity and for creation of frustrated lattices to explore the 
possibility of unconventional quantum phases p} ■ It may 
also be possible to achieve reliable prototypes for quan- 
tum computing and information processing [8j]. Recent 
experiments involving dipolar atoms [9( are expected to 
lead to the realization of new phases of matter such as 
the supersolid (SS) 0, 111! , whose existence is still not 
unambiguous in solid 4 He [12j. 

The MI-SF transition in lattice bosons is controlled 
by the competition between the interaction among them 
and their kinetic energy [l3j]. However, the translational 
symmetry breaking confining potential, that is present 
in experiments, would lead to inhomogeneous phases 
and even phase coexistence for appropriate parameters. 
There have been various proposals to observe such struc- 
tures [14|. A clear indication of it will be their excita- 
tion spectra, which will be different from their homoge- 
neous counterparts and the way the system behaves un- 
der time-dependent perturbations of the trap [la ]. Fur- 
ther, the study of the collective modes has become an 
important experimental tool to analyze correlation prop- 
erties of ultracold quantum gases in a trap 0, [l6| • With 
this in mind, we systematically explore the excitations of 
parabolically trapped lattice bosons from deep SF phase 
to insulating phase. Our method also provides a detailed 
microscopic understanding of the behavior of the trapped 



SF under perturbations which, hitherto, has been studied 
using the Gross-Pitaevskii equation (GPE) or equivalent 
hydrodynamic approach [ItT ]. These methods are valid 
only when the SF fraction is very large; they fail to cap- 
ture the modified physics both in the homogeneous SF 
near the MI boundary and in the inhomogeneous phase 
coexistence regime. 

The harmonically trapped alkali atoms with short 
range interactions in an optical lattice can be modeled 
by the effective Bose-Hubbard model (BHM) P,[H: 
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The three terms represent the lattice kinetic energy, the 
on-site interaction, and the harmonic trap potential re- 
spectively. The operator aj creates a boson at site i, 
hi = a\hi is the boson number operator, is the dis- 
tance of the site i from the minimum of the trap poten- 
tial with frequency to; t,U are the nearest-neighbor hop- 
ping and the on-site Hubbard interaction respectively, m 
is the mass of the atoms, and S represents the nearest 
neighbors of the site i. This Hamiltonian without any 
confining potential has been studied extensively using a 
variety of methods such as the mean field theory [l3j . 
quantum Monte-Carlo [3], perturbation series [h| and 
density matrix renormalization group (DMRG) [2CJ, |21 1 . 
In the absence of trap, the ground state is an incompress- 
ible MI with an integer number of atoms at every lattice 
site when the interactions dominate and the particle den- 
sity is commensurate with the lattice. When the kinetic 
energy dominates, and in general when the density is 
incommensurate, a SF ground state is obtained. The 
low-lying excitations in the MI state are gapped particle- 
hole excitations. The SF phase has gapless, acoustic, 
Bogoliubov quasiparticles. For large enough t/U, the 
excitations can be described by the discrete non-linear 
Schrodinger(DNLS) equation (the lattice version of the 
GPE) 17], which becomes progressively inaccurate as 



2 



one moves towards the MI-SF boundary In pres- 

ence of a trap one encounters an inhomogeneous SF phase 
and at smaller values of t/U, the renormalized chemical 
potential gives rise to alternating shells of SF and MI re- 
gions. In absence of numerically exact calculations, which 
are not hindered by limited system size or small num- 
ber of bosons, a (linearized) time dependent variational 
mean field analysis is an appropriate tool to explore the 
excitations of these systems; it gives quantitatively cor- 
rect results [1S|] in homogeneous case. Hence, we use 
the variational Gutzwiller mean field approach to study 
the ground state and use a modified time-dependent 
Gutzwiller ansatz to study the low-lying excitation en- 
ergies. In particular, we concentrate on the excitations 
in the homogeneous phase near the MI-SF boundary, the 
low-lying spectra in the presence of a static trap as a 
function of parameters, and the dynamics of the system 
in response to time-dependent perturbations of trap po- 
tential. 

We calculate the inhomogeneous ground states of the 
above model, in general in d dimensions, for a given 
set of parameters and the chemical potential /i by mini- 



mizing 



H — /j,N 



^ ) with a Gutzwiller wavefunction 



\ty) = Yli J2 n /™ I 71 ' *) with respect to the variational pa- 
rameters fn \ where \n, i) is the Fock state with n par- 
ticles at site i and N = is the total particle num- 
ber operator. The excitation energies above the ground 
state are, then, obtained from the real-space dynami- 
cal Gutzwiller approach with the variational parameters 

(i) 

fn _ being time dependent that was introduced in Ref. 
22] and modified for the calculation of excitation spec- 
tra in Ref. [111. Minimization of the effective action 
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Here fa — (\E , |cij|\E r ) is the local condensate (SF) or- 
der parameter and r is time. The oscillation frequencies 

of the small amplitude fluctuations Sfjfi (r) around the 

yd) 

ground state /A give the excitation spectrum. Normal- 
ization of the wavefunctions, ^ n \fn\T)\ 2 = 1, at each 
site is enforced using a Lagrange multiplier A^. 

Firstly, we consider the excitations in the homogeneous 
system without a trap. For the MI phase with n a bosons 
per site, the ground state is defined by f„ = 5 n ,n - 
The resulting linearized equations lead to particle(p) and 



hole(h) excitations with dispersions 
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where ± corresponds to p(h) excitations and ek 
— 2t 5Zj=i cos(fej). This is in agreement with the re- 
sults obtained from the slave-boson mean-field theory 



and random- phase approximation 24 1 



As one ap- 
proaches the phase boundary from the MI side, the en- 
ergy gap for p(h) excitations gradually decreases, van- 
ishing at the transition. The phase boundary at which 
the MI-SF transition occurs is obtained from the condi- 
tion of vanishing energy-gap, e p (^)(k = 0) = 0, and is in 
agreement with the calculations using second order per- 
turbation theory in t/U [19]. Deep in the SF phase (i.e., 
when t > £/), and for large SF fraction (\fa\ 2 ~ ni), the 
wavefunction at each site can be represented by coherent 
states, i.e., fn = fa/e'^^ I 2 / \fn\ and Eq.(2) reduces to 
a DNLS equation for the classical field fa: 
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In the absence of a trap, this gives gapless acoustic mode 
with sound velocity c s = \ fa\ay2tU /h, where a is the lat- 
tice spacing. The MI-SF transition near the Mott lobe 
at a commensurate density no can be understood via a 
simple three-state variational ansatz at site i: \ipi) = 
fi | no — l)+/o |^o) + /2 \no + !)■ This captures the build- 
up of number fluctuations [IH, and hence the phase co- 
herence as the system makes a transition to the SF phase 
at a critical coupling (U/2td) c = (y/no + V n o + 1) ■ In- 
troducing time-dependent fluctuations in the variational 
parameters lead to the Bogoliubov spectrum in the SF 
phase with the sound velocity c s given by, 
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where a = (y/n~o + yfno + 1)" = U/(2tdcos29). Interest- 
ingly, at the boundary where 9 = 0, the sound velocity is 
finite, although the SF order parameter <j> vanishes. This 
was first pointed out in Ref. [26( with no ^ 1 and an ef- 
fective relativistic GPE has been proposed to capture the 
dynamics close to the MI phase. For small momenta, the 
(gapped) amplitude mode [26[ in the SF has the disper- 
sion w(k) = VA 2 + c 2 k 2 with A = td sin 29y/a 2 - 1 and 
c 2 = (t 2 d/8) [9 + (4a 2 - 13) cos 2 29] in d-dimensions. 
At the boundary, this mode becomes degenerate with 
the gapless, Bogoliubov mode. However, deep in the SF 
phase three state variational ansatz is insufficient and 
we increase the number of Fock states at each site and 
numerically solve the linearized equations for the fluc- 
tuations Sfn^ to obtain the collective mode with mo- 
mentum k. The dispersion of the lowest acoustic mode 
in the SF phase is shown in Fig.l for different coupling 
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strengths. We compare our result with those obtained 
from DNLS and find good agreement deep in the SF 
regime (\4>i\ 2 ~ n,), whereas DNLS fails near the MI 
boundary (see Fig.l). 

Main advantage of our method is that it is simple 
and efficient enough to implement in confined geometries 
and for inhomogeneous systems. Confining traps which 
are present in experiments can be approximated by a 
harmonic oscillator potential V(rj) = ^mu> 2 r 2 . Mini- 
mization of free energy gives, in general, inhomogeneous 
ground state that show coexisting SF and MI phases 
(with step-like structure) for certain parameters. Com- 
pressible edge forms on the boundary of such droplet. 
Having found the ground state for a given /z, we cal- 
culate the low-lying spectrum by solving the system of 
linearized equations (Eq.(2)). Recently collective excita- 
tions of BHM in a ID harmonic trap have been calculated 
by exact diagonalization of the Hamiltonian for small sys- 
tem size However, the complexity of the problem 
grows exponentially with increasing number of the lat- 
tice sites and boson filling fractions. Numerically efficient 
T-DMRG methods dUEl 

are restricted to ID systems 
with short range interactions, whereas our method can 
be applied to higher dimensions, where the mean field 
method is expected to work better. In Fig. 2 we show 




FIG. 1: Excitation spectrum of the uniform BHM with 
fx/U = 0.5, a) for t/U = 2 and b) for t/U = 0.087. For com- 
parison dotted lines represent excitation spectrum obtained 
from DNLS. Momentum k is in units of n/a. 

a few low-lying excitations of bosons in a ID harmonic 
trap, as a function of t/U for a fixed fi/U. We clearly 
notice three different regimes by changing the interaction 
strength t/U. In the deep SF region (i.e., t/U 1), exci- 
tations are Bogoliubov quasiparticles and their energies 
asymptotically match with those obtained from DNLS. 
In this regime the hopping parameter t is equal to the 
kinetic energy h 2 /(2ma 2 ). For large t/U, the nth excita- 
tion branch approaches to y/ n(n + l)ta 2 m asymptotically 
and is in agreement with the hydrodynamic modes of 
a dilute Bose gas in ID harmonic trap without a lat- 
tice [28| . As the system enters the correlated regime (for 
smaller t/U), which is close to the phase boundary, the 
excitation energies start to deviate from those obtained 
from DNLS and many avoided level crossings occur in 
the energy spectrum. We notice the clear signature of 
particle and hole type excitations deep in the MI phase 
(t/U -C 1). In contrast to the homogeneous BHM, where 
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FIG. 2: Few low- lying excitations of the BHM in a trap with 
mu> 2 a 2 /U = 0.16 and [i/U = 0.6 (solid lines). For comparison 
the dotted lines denote the energies obtained from DNLS for 
the same parameters. 



p(h) energies are site-degenerate, in case of harmonically 
trapped BHM, they split and become site-dependent in 
the insulating phase. It is easy to understand these p(h) 
excitation energies in the atomic limit of BHM. For t = 
and n < U the system forms a droplet of size 2n max 
with all sites i for — n max < i < n max filled with one 
particle per site, where n max being the integer part of 
\/2/i/(mw 2 a 2 ). Thus in presence of the trap, the energy 
for adding or removing a particle, in general, depends on 
the site index i, and is given by E p = U — \i + ^mui 2 a 2 i 2 
and Eh = \i — \muj 2 a 2 i 2 . At the edge of the trap (at 
i = ±n max ), E p « U and Eh ~ 0. Similarly, it costs 
very little energy to create a particle at the empty site 
just outside the droplet (i.e., i = ±(n max + 1)). This 
leads to the formation of the gapless edge at finite t. 
In SF phase, among the low-lying collective excitations 
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FIG. 3: Time evolution of the breathing mode Ax 2 in a trap 
with muj 2 a 2 = 0.5, p = 0.8, a) for t/U = 0.23, and b) for 
t/U — 0.1. Length and time are measured in units of a and 
h/U respectively. 

of harmonically trapped lattice bosons, dipole oscillation 
describing the center of mass motion has the lowest en- 
ergy. According to Kohn's theorem, dipole excitation 
energy in a harmonic trap is ftw that is independent of 
two-body interaction [13]. In case of BHM the dipole 
excitation energy deviates from this universal value due 
to the presence of the optical lattice. For large system 
size and deep in the SF phase we find that it approaches 
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\j2brruoa, in agreement with Kohn's theorem. In exper- 
iments, dipole mode can be excited by slightly shifting 
the atomic cloud from the center of the trap by imposing 
a perturbation of the form V per t ~ Xihi. 

Having studied the spectra in the presence of a static 
trap, we now concentrate on the effect of time-dependent 
perturbations of the trap on the dynamics of BHM by 
solving the full time-dependent Gutzwiller equation. As 
an example, we focus on an experimentally relevant col- 
lective mode of ID confined system, the breathing mode, 
where the center of mass of the atomic cloud remains 
fixed but its size oscillates. For harmonically trapped di- 
lute Bose gas in quasi ID regime, the frequency of this 
oscillations is v^w [28| . Experimental study of this mode 
has become an important tool for the investigation of the 
correlation properties of quasi ID quantum gases, partic- 
ularly when ID Bose gas enters the strongly correlated 
fcrmionized regime [28]. We study the breathing mode 
as a linear response to a perturbing potential of the form 
xjfii. Oscillation frequencies of both dipole 
and breathing modes obtained from the full dynamics 
agree with those obtained from linearized time depen- 
dent method. In Fig. 3, time evolutions of the collective 
coordinate Ax 2 = [(xf (r)) — (x 2 (0))] are shown for 
two different values of t/U , under the same perturbation 
strength. We observe that the response of the system 
becomes smaller as it approaches the insulating regime. 
The MI phase hardly responds to perturbations due to 
the energy gap in the spectrum [29(. In the SF phase 
with small t/U (close to the MI boundary), the breath- 
ing mode oscillations show beats instead of a single fre- 
quency. As noted earlier, in this regime excitation ener- 
gies are very close to each other and many avoided level 
crossings occur (see Fig. 2). Due to this reason mode 
coupling takes place leading to beats in the dynamics of 
the collective modes. 

In conclusion, we develop a new time-dependent vari- 
ational method to obtain the spectra of correlated lat- 
tice bosons in an (arbitrary) confining potential. In the 
homogeneous phase, we find that, at commensurate fill- 
ing and at the SF-MI boundary, the sound velocity does 
not vanish even though the SF order parameter van- 
ishes. Low-lying excitations of BHM in a harmonic trap 
are obtained, which validates the GPE approach asymp- 
totically, while deviating significantly in the correlated 
regime. The dynamics of the collective coordinates under 
time-dependent perturbations reveals decrease in their 
amplitudes as well as appearance of beats in the strongly 
correlated regime. These results are relevant to the ex- 
perimental observation of coexisting MI-SF phases and 
measurements of their dynamical properties. 
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